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Abstract. We study the Ho- emission from jets using two-dimensional axisymmetrical simulations. We compare the emission 
obtained from hydrodynamic (HD) simulations with that obtained from magnetohydrodynamics (MHD) simulations. The mag¬ 
netic field is supposed to be present in the jet only, and with a toroidal configuration. The simulations have time-dependent 
ejection velocities and different intensities for the initial magnetic field. The results show an increase in the Ha emission along 
the jet for the magnetized cases with respect to the HD case. The increase in the emission is due to a better collimation of 
the jet in the MHD case, and to a small increase in the shock velocity. These results could have important implications for the 
interpretation of the observations of jets from young stellar objects. 
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1. Introduction 

Collimated outflows are observed in a variety of astrophysical 
objects, with typical spatial scales ranging from ~ 1 pc for jets 
from young stellar objects (YSOs) up to several megaparsecs 
for extragalactic jets. All of these jets seem to be associated 
with accretion disks, which suggests the existence of a scale- 
independent physical mechanism responsible for the ejection 
and collimation of these outflows. The presently more accepted 
models are the magnetocentrifugal models (Blandford & Payne 
fT^ Uchida & Shibata [T985l . in which the ejection is driven 
by the presence of a dynamically important magnetic field in 
the accretion disk-central object system, and the collimation of 
the jet is due to the toroidal component of the magnetic field, 
which is able to collimate the outflows by pinching forces. The 
toroidal magnetic field is generated by the twisting of the mag¬ 
netic field due to the rotation of the system. The region where 
this process acts is too close to the central object to be resolved 
observationally, and one possible way to obtain some insight 
into this region is by studying the properties of the outflows. 

In particular, a lot of progress has recently been made re¬ 
garding observations of the outflows from YSOs (see the re¬ 
view by Reipurth & Bally 1200 It . The caracteristic spectral 
emission of these objects is believed to come from the region 
behind the shock, from recombination of the ionized gas (for 
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the hydrogen lines) and electron excitation (and de-excitation) 
within ions (Schwartz ll975t . and the typical knot structure visi¬ 
ble along the jet could be interpreted as due to a time-dependent 
ejection from the young stellar object (YSO) (Reir)urth ll989l 
Raga et al|T990j. 

In the past few years, several authors have studied the ef¬ 
fect of the magnetic field on the dynamical evolution of HH 
objects (e. g. O’Sullivan & Rav I200()l Stone & Hardee EOOOl 
Cerqueira et al. I1997> . Frank et al. (I1999t showed that am- 
bipolar diffusion could be important to smear out the magnetic 
field ejected with the jets, but only on timescales comparable 
or larger than the dynamical timescales of the jets. The deter¬ 
mination of a magnetic field in the outflows would represent 
a test for the magnetocentrifugal mechanism. In particular, the 
presence of a dynamically important toroidal component of the 
magnetic field would represent an indirect proof of this mech¬ 
anism. However, direct observations of magnetic fields in jets 
are very difficult, and there is yet no clear observational deter¬ 
mination of the magnetic field intensity in jets. 

Several hydrodynamic (HD) simulations with calculations 
of the spectral emission of HH objects have been presented in 
the past (e. g. Blondin et al. 119901 Raga ll99?b . but such cal¬ 
culations have never been published for the magnetized case. 
The simulations with magnetic fields have concentrated on the 
dynamical aspects and the evolution of the jet rather than on 
obtaining predictions of the emitted spectrum. These simula¬ 
tions usually include a radiative cooling rate (given, e.g., by 
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the coronal cooling function of Dalgarno & McCray wn\ . 
and different magnetic field configurations. The main features 
found in magnetized jets with respect to HD jets are the pres¬ 
ence of a “nose cone”, better collimation, an increase in the 
density along the direction of propagation, and some effects 
on Kelvin-Helmoltz and Rayleigh-Taylor instabilities (result¬ 
ing in changes in the leading bow shock, e. g. Todo el al ll993l 
Cerqueira & de Gouveia Dal Pino 1 199^ . 

With respect to the emitted spectrum, Hartigan, Morse & 
Raymond J1994t compared observed and predicted emission 
lines ratios (using plane-parallel shock models) to find an upper 
limit of 30 for the magnetic field of the jet. Cerqueira & 
de Gouveia Dal Pino J2001a> . using a semiempirical formula 
(valid for a shock velocity between 20 and 80 km s“') found 
the ratio between the Ha emission of a magnetized and a non- 
magnetized jet. They obtained that the Ha emission increases 
due to the presence of a magnetic field, and that the dominant 
cause of this increase is the toroidal component of the magnetic 
field. 

Therefore, no direct calculation of the effect of the mag¬ 
netic field on the emission from a jet has yet been made. Trying 
to fill this gap, we have carried out 2D, axisymmetric, MHD 
simulations of jets, looking for the differences in the predicted 
Ha emission from magnetized jets with respect to the hydro- 
dynamic case. The paper is organized as follows. In section 
2, we explain in some detail the numerical algorithm that has 
been used, the initial conditions of the simulations and the ap¬ 
proximations used to calculate the emission. In section 3 we 
summarize and discuss the results obtained, and in section 4 
we draw our conclusions. 


2. Simulations 

2.1. Numerical method 


To study the Ha emission from YSO jets, we have carried out 
a set of 2D, axysimmetrical simulations using a MHD code, 
modified to include the cooling and the evolution of the hy¬ 
drogen ionization fraction. We solve the following set of equa¬ 
tions: 


^ -H V ■ (pv) = 0 
ot 


( 1 ) 


^ -I- V ■ (pvv + p,o,I - BB) = 0 
ot 


( 2 ) 


de 

-H V ■ ((e H- pto,)v ■ 
at 


(V B)B) = 


-n^A 


( 3 ) 


^ -I- V ■ (vB - Bv) = 0 
at 


( 4 ) 


dufjo 

dt 


-H V ■ (tiffov) - tiH+ngaiT) - riffarieCiT) 


( 5 ) 


where p is the mass density, v is the velocity vector, ptot - 
Pgas + B^/2 is the (magnetic + thermal) total pressure, I is the 
identity matrix, B is the magnetic field normalized with respect 
to V4^, e is the total energy defined as e = pgas + \pv^ + 5 ^^ 
(with y - 5/3), and finally A is the cooling function. These 


equations represent the conservation of mass Q, momentum 
(13, energy 0 and magnetic flux 0. Eq. 0 represents the 
evolution of the hydrogen neutral fraction. This equation is 
coupled to the others by the cooling function A, which is a 
function of the ionization fraction. In Eq. 0 : 

a{T) = 2.55 ■ 10^'^(7’/10'^)“*’™ 
ar) = 5.38 ■ 10^"e-^/Vf 


are the recombination and collisional ionization coefficients 
(Cox ll97(?l . and T is the gas temperature. 

To solve this equation system, we use a second order up¬ 
wind scheme, which integrates the MHD equations using a 
Godunov method with a Riemann solver. The Riemann prob¬ 
lem is solved using primitive variables and the magnetic held 
divergence is maintained close to zero using the CT method 
(Toth I2000t . The algorithm is similar to the one of Falle, 
Komissarov and Joarder J1998t . except that to preserve a small 
divergence of the magnetic held we use a constrained trans¬ 
port method. However, the magnetic held divergence is au¬ 
tomatically equal to zero in all our axysimmetrical simula¬ 
tions because we use a toroidal magnetic held configuration. 
To treat correctly the source terms in cylindrical coordinates, 
these terms are averaged on the cell volume (see the discussion 
in the appendix of Falle ll99Tl . The code was tested with one- 
and two-dimensional tests (Ryu et al. 119951 Toth KOOOl also see 
De Colle l200^ ). 

To solve Eq. 0 , we first integrate the equation without the 
cooling term finding a new value for the energy e', and after¬ 
wards we use this new energy value to integrate the equation 

de' 0 j T 

— = -n^A(T) = -n^Ao —, (6) 

dt To 


where n is the number density (assumed to be constant within 
the timestep) and we approximate A(T) in a locally linear form. 
This equation has an exact solution for the gas pressure: 


Pgas — P0,gas£ 


—2AoW// 


( 7 ) 


where hh is the total (ionized more neutral) hydrogen density, 
T is the gas temperature, ks is the Boltzmann constant and At 
is the time step. A floor of 1000 K is used as the minimal tem¬ 
perature value. 

Tests for adiabatic and radiative magnetized and HD jet 
propagation (reproducing the results of Frank et al. IT^ are 
reported in a previous paper (De Colle & Raga l20d?i . 


2.2. Cooling term 

We have used a non-equilibrium cooling function considering 
the energy loss due to collisional excitation of oxygen, radiative 
recombination of hydrogen, collisional ionization of the hydro¬ 
gen and excitation of Lyman-alpha line (Biro et al. ll995t . If one 
computes an coronal equilibrium cooling curve using this sim¬ 
plified scheme, one obtains a cooling function which is similar 

' A detailed description of the free distributed numerical code with 
several different one, two and three dimensional numerical tests will 
he soon available at www.astroscu.unam.mx/~fdecolle 
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Fig. 1. Numerical density for the models HD (hydrodynamic), weak MHD (fi — V), and strong MHD (fi - 0.4), after a 500 yrs 
integration time. The right bar gives the logarithm of the numerical density (in cm“^) The axes are labeles in pixels, and the 
displayed domain has a physical size of {L~,Lr) = (3 x 10*^, 3 x lO'®) cm. A pixel correspond to a physical size of 1.66 x lO'"^ cm. 


to the one of Dalgarno & McCray (I1972t . which is usually used 
in MHD jets simulations. For temperatures between 10"^ K and 
10^ K it is necessary to use a non-equilibrium function instead 
of the simpler equilibrium cooling function. However, using 
our non-equilibrium cooling function, we find that the dynam¬ 
ical evolution of the jet and the main features (such as better 
collimation, presence of a small nose cone, etc.) are very simi¬ 
lar to the ones of previous simulations (obtained with a coronal 
equilibrium cooling function). 

To obtain the Her emission we consider the contributions 
from the radiative recombination cascade (Aller|l984j and col- 
lisional excitations from the n = 1 state (Giovanardi and Palla 
119891 . The Ha emission is then obtained as a function of the 
local temperature of the gas T, and of the neutral and total den¬ 
sities riH- 


2.3. Initial conditions 

In our numerical simulations, we use a uniform grid of ax¬ 
ial/radial size m,Lr)-{3 X lO'^, 3 x lO'®) cm. The resolution is 
of 1.67 X 10'"^ cm, corresponding to 1800 x 180 cells along the 
Z- and r-axes, respectively. 


In the models, the atomic jet has an initial radius of 6 x 
lO'^ cm (corresponding to 30 cells) and is moving in the z- 
direction with a mean velocity vo = 200 km s“' (on this mean 
velocity, we superimpose a sinusoidal velocity variability, see 
below). The initial jet density in all cases is tije, = 500 cm“^, 
the ambient density is uniform, with a value namh = 100 cm~^, 
and the ambient sound speed is Cj = 10 km s * (corresponding 
to r = 10000 K). The edge of the jet and the ambient medium 
are initially in radial pressure equilibrium. The initial jet tem¬ 
perature has a radial profile which is discussed below. 

The jet evolution is followed for a maximum of 500 yr as 
it propagates along the z-direction. We use reflection boundary 
conditions at z=0 and r=0, and open boundary conditions for 
the outer boundaries of the x and z-axes. 

To obtain a number of knots along the direction of prop¬ 
agation of the jet, we impose a sinusoidally variable ejection 
velocity of the form: 


v(f) = Vo 


/ 27 rf\ 

II H- A sin — I 


( 8 ) 


where A = 0.25 and t = 50 yr are the amplitude and the period 
of the perturbation in the velocity. 
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Fig. 2. Axial cuts of the density (top), hydrogen ionization fraction, axial velocity and temperature (bottom) for the different 
models. The solid lines correspond to the HD model, the dashed lines correspond to the weak magnetic field (j3 - 1) model, and 
the dotted lines correspond to the strong magnetic field (/? = 0.4) model. 


Following Lind et al. dnn, we use the simple profile for 
the toroidal magnetic field: 


B(r) = 


B 




0 < r < R„ 


Bm-r r,n<r <Ri 


( 9 ) 


0 


Rj < r 


To ensure initial hydromagnetic equilibrium in the radial di¬ 
rection, it is necessary that the magnetic and thermal pressures 
satisfy the radial equilibrium equation: 


dp ^ B d(rB) ^ 
dr r dr 


( 10 ) 


Integrating this equation, using the magnetic field profile given 
by Eq. it is possible to find the thermal pressure profile 
(Lind et al. fTOS^ O’Sullivan & Ray|2000j: 


p{r) 


"I" d2 y)Pamb 0 < r < Rff 


I m — r < R j 
Rj<r 


( 11 ) 


^Pamb ^tr. 

Pamb 

where a is the integration constant of Eq. CoJ, which is con¬ 
nected to Rm by the relation; 

\2 


a - I 


1 

(dm 


( 12 ) 


This is the magnetic field prohle used by Lind et al. ifT^ . 
O’Sullivan & Ray | |2000| | and Erank et al. (I1998t . and is similar 
to the one Stone & Hardee J2000I who use a slightly different 
form for the magnetic held prohle). 

The free parameters in these formulae are the values of R,„ 
(or equivalently a). We choose R,„-Q.6Rj. In this way the a 
parameter (obtained from Eq. d, is always positive and the 
thermal pressure (obtained from Eq. CD is also positive. 

We choose a toroidal held geometry for two reasons. The 
hrst one is that a toroidal magnetic held is necessary to colli¬ 
mate the jet during the ejection-collimation process, and would 
then also have to be present in the propagating jet. The sec¬ 
ond reason is that in shocks, the “important” component of the 
magnetic held is the one parallel to the surface of the shock 
(which is the toroidal component in our geometry). In fact, all 
of the previous simulations conhrm this idea. Stone & Hardee 
and Erank et al. ifT^ obtained relevant differences in 
MHD jets (with respect to the HD case) with a strong toroidal 
component of the magnetic held, and less important differences 
for a longitudinal magnetic held (see also Gardiner et al. l2000t . 
Cerqueira & de Gouveia Dal Pino J2001all2001b> also showed 
that larger effects in the calculated Ha emission are also ob¬ 
tained for toroidal held conhgurations. 
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0 MO 1000 1500 

Fig. 3. Ha emission maps for the three computed models. The frames correspond to the HD, weak MHD ij3 - 1), and strong 
MHD ijS - 0.4) models (up to bottom), after 500 yrs integration time. The right bar gives the logarithm of the Ha emission, 
normalized to the maximum value of each frame. The axes are labeles in pixels, and the displayed domain has a physical size of 
{L-,Lr) = (3 X lO'^, 3 X lO'®) cm. A pixel correspond to a physical size of 1.66 x 10'"^ cm. 


We have computed models with different magnetic field in¬ 
tensities, ranging from - 0 to - 93.2 fiG (see equation 
13, corresponding to p - oa to p - 0.4, (where p - IPgasIB^). 
In all of our simulations, the ambient medium is not magne¬ 
tized. The sound velocity and the Alfven velocities, at the cen¬ 
tre of the initial cross section, ranges from Cj_o -Cs,ambl VS a; 
4.5 km/s and Ca = 0 (for the p - oo case) to Cj_o ~ 10 km/s and 
ca ~ 17.3 km/s (for the p - 0.4 case). 

3. Results 

3.1. General Results 

The stratification of the numerical density for simulations with 
different values of the magnetic field intensity are shown in Fig. 
1 for a f = 500 yr time integration. The top frame shows the HD 
simulation, the central frame a weak magnetic field simulation 
(with P=l, B„, = 58.9 fiG) and the bottom frame a strong mag¬ 
netic field simulation (P-0.4, = 93.2 /uG). 

The results show an increasing collimation for increasing 
toroidal magnetic fields. The knots are due to the imposed vari¬ 
able injection velocity (see Eq.|3- The head of the jet presents 


a very complicated structure, due to the presence of instabil¬ 
ities (see, e. g., Blondin et al. 1990) and to the interaction of 
successive knots which catch up with the head. 

An important result is that the Ha emission of the leading 
head is strongly dependent on p. This is a result of the fact that 
our magnetized jet simulations develop ’’nose cones”. These 
structures are of somewhat dubious reality, as they might dis¬ 
appear in 3D jets without perfect axisymmetry (see Cerqueira 
& de Gouveia Dal Pino l200rbt . 

The structure of the head is strongly time-dependent, and 
will not be considered in the following analysis. In Fig.|2lwe 
show the density, the ionization fraction, the axial flow velocity 
and the temperature as a function of position along the symme¬ 
try axis for the three computed models. This figure does not 
include the region of the head of the jet. 

Some effects due to the presence of the magnetic field are 
evident. We see that the on-axis (r = 0) density increases ~ one 
order of magnitude from the HD to the P-Q.A model. The ini¬ 
tial temperature is bigger for the magnetized cases due to the 
different pressure profiles present in the initial configuration of 
the jet (which correspond to temperature profiles, as the den- 
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sity is assumed to be constant, see section 2). However, within 
the internal working surfaces we obtain lower temperatures as 
a function of increasing magnetic held intensity (similar results 
were obtained by Stone & Hardee I2000> . The ionization frac¬ 
tion in the central regions of the knots decreases as a function 
of increasing magnetic held strength (which is consistent with 
the temperature decrease described above). 

The separation between the two shocks of the knots in¬ 
creases both as a function of distance from the source and of 
increasing magnetic held intensity. This latter effect is due to 
the fact that the magnetic held tends to stop the lateral expan¬ 
sion of the material within the knots, so that more material re¬ 
mains in the region close to the symmetry axis. With respect to 
the velocity, we see that the shock velocities (which correspond 
to the velocity jumps observed at the edges of the knots in Fig. 
2) increase as a function of increasing magnetic held. Finally, 
the propagation velocity of the head of the jet increases sub¬ 
stantially as a function of increasing magnetic held (as can be 
seen in Fig. 1). This effect is also due to the increased collima- 
tion obtained for the MHD jets. 

3.2. Intensity maps 

From our models, we have computed Ha emission line maps 
from the models. These maps are obtained integrating the emis¬ 
sion coefficient along lines of sight, and are computed assum¬ 
ing an optically thin jet moving parallel to the plane of the sky. 
In Fig. 3 we show the Ha emission line maps predicted from 
the models presented in Fig. 1. It is evident that the Ha emis¬ 
sion is more concentrated to the jet axis in the magnetic jet 
models than in the HD jet. 

To have a quantitative value of the Ha emission from the 
knots, we integrate the emission over regions around each knot. 
To this effect, we integrate the emission within circular di¬ 
aphragms with a radius of 7 x 10^^ cm centered on the succes¬ 
sive knots. In fact, the regions outside these diaphragms have 
negligible emission (see Fig. 3). 

In Fig. 4, we show the results obtained for the different 
models. The points represent the ratio of the Ha luminosi¬ 
ties of the successive knots with respect to the Ha lumninos- 
ity of the corresponding knot in the HD simulation. For the 
hrst knot (located at z ~ 4-Rj) we hnd that all of the magne¬ 
tized jet models have almost the same Ha luminosity as the 
HD jet model (i. e., LnaiMHD)!LnaiHD) » 1). The value of 
LnaiMHD)ILnaiHD) increases as a function of z in the second 
and third knots (located at z 8 and l4Rj) and then decreases 
again for larger z values (i. e., for the fourth and hfth knots, 
located at z ~ 19 and 25Rj). 

From this Fig. 4, we also see that a model with P - 5 
(B - 26.3jL(G) produces knots with Ha luminosities which are 
basically indistinguishable from the HD simulation. For higher 
magnetic helds (i. e., for lower values of P), knots with in¬ 
creasingly larger Ha luminosities (compared to the HD case) 
are obtained. The largest increase is found for the knot at po¬ 
sition z ~ I4rj (» 8.5 X 10'® cm) These results agree at least 
qualitatively with the ones shown in Fig. 3 of Cerqueira and de 
Gouveia Dal Pino J2001al l. who hnd increases in the Ha lumi- 



Fig.4. Ratios of the Ha luminosities of the knots along magne¬ 
tized jet models with respect to the Ha luminosity of the corre¬ 
sponding knots in a HD jet, as a function of knot position along 
the jet axis. These luminosities have been obtained by integrat¬ 
ing the Ha intensity over circular diaphragms centered on the 
successive knots (see the text and Fig. 3). 

nosity by factors of ~ 3-4 with respect to the HD case for a 
P - \ magnetized jet model. 

4. Conclusions 

We have presented 2D numerical simulations of HD and mag¬ 
netized, variable jets propagating in a homogeneous medium. 
Models with increasing magnetic helds show an increase in the 
Ha luminosities of the successive knots (which correspond to 
internal working surfaces which result from the injection ve¬ 
locity variability). This result conhrms the work of Cerqueira 
& de Gouveia Dal Pino (I2001all2001bt . who estimated the Ha 
luminosity of the clumps along MHD jet simulations using a ht 
to predictions of plane-parallel shock models. 

Somewhat surprisingly, our work presents the hrst pre¬ 
dictions of emission line maps from MHD HH jet models. 
Therefore, our calculations for the hrst time show the emission 
line morphologies that would be expected for such models. 

We hnd that the Ha emission of the leading head of the jet 
differs quite strongly between the HD and MHD cases. This 
is a result of the fact that our simulations develop extended 
“nose cones” (of somewhat dubious reality, as these structures 
might disappear in 3D jets without perfect axisymmetry, see 
Cerqueira & de Gouveia Dal Pino 1200Ibt . 

For the knots along the jet, we hnd that for increasing mag¬ 
netic held strengths we obtain emission structures with stronger 
peaks towards the symmetry axis. This can be clearly seen in 
Fig. 3, in which the p - 0.4 model has knots which are domi¬ 
nated by an elongated emission component along the jet axis. 

This different type of knot morphology is interesting in 
terms of observations of HH jets. It has been a long-standing 
fact that while some HH jets (notably HH 111, see e. g. 
Reipurth et al. I1996> show compact knots with “bow shock¬ 
like” morphologies which resemble the predictions from HD 
variable jet models (see Masciadri et al. I2002> . other HH jets 
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(e. g., HH 30, see Lopez et al. ll995l l have emission knots with 
axially elongated structures. This second kind of morphology 
could not be modeled successfully in terms of variable HD jet 
models, and suggested the presence of a different mechanism 
for knot formation. We now find that variable jets with a strong 
enough toroidal magnetic field do lead to the formation of axi¬ 
ally elongated knots, which in principle might be used to model 
objects such as HH 30. 

The present paper is limited to a study of the effect of a 
toroidal magnetic held on the Ha emission of variable jet Hows. 
In a future paper, we will present a study of a more extended set 
of emission lines. From such a study, we will attempt to pro¬ 
duce a set of line diagnostics which could be used to estimate 
the magnetic Held strength along observed HH jets. 
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